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Nickel-Rich Outflows Produced by the Accretion-Induced 
Collapse of White Dwarfs: Lightcurves and Spectra 
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ABSTRACT 

The accretion-induced collapse (AIC) of a white dwarf to form a neutron star can 
leave behind a rotationally supported disk with mass of up to ~ 0.1 M©. The disk is 
initially composed of free nucleons but as it accretes and spreads to larger radii, the free 
nucleons recombine to form helium, releasing sufficient energy to unbind the remaining 
disk. Most of the ejected mass fuses to form 56 Ni and other iron group elements. We 
present spherically symmetric radiative transfer calculations of the transient powered 
by the radioactive heating of this ejecta. We estimate the ejecta composition using 
nucleosynthesis calculations in the literature and explore the sensitivity of our results 
to uncertainties in the ejecta kinematics. For an ejecta mass of 10 -2 M© (3 x 10~ 3 M©), 
the lightcurve peaks after < 1 day with a peak bolometric luminosity ~ 2 x 10 41 ergs s _1 
(~5x 10 40 ergss _1 ); the decay time is ~ 4 (2) days. Overall, the spectra redden with 
time reaching U — V ~ 4 after ~ 1 day; the optical colors (B — V) are, however, 
somewhat blue. Near the peak in the lightcurve, the spectra are dominated by Doppler 
broadened Nickel features, with no distinct spectral lines present. At ~ 3 — 5 days, 
strong Calcium lines are present in the infrared, although the Calcium mass fraction is 
only ~ 10~ 4 - 5 . If rotationally supported disks are a common byproduct of AIC, current 
and upcoming transient surveys such as the Palomar Transient Factory should detect 
a few AIC per year for an AIC rate of ~ 10~ 2 of the Type la rate. We discuss ways 
of distinguishing AIC from other rapid, faint transients, including .la's and the ejecta 
from binary neutron star mergers. 

Key words: supernovae - nucleosynthesis, abundances - accretion disks 



1 INTRODUCTION 

As an accreting white dwarf (WD) approaches the Chan- 
drasekhar mass, it suffers one of two fates: (1) collapse 
to form a neutron star (NS) or (2) explosion as a Type 
la supernova (SN). Which of these end states is realized 
in dif ferent astrophysical scen arios is not fully understood 
(see ICanal fc Gutierrez! 1 19971 for a review). The composi- 
tion of the WD is, however, a critical factor: the relatively 
low densities required for electron captures on Ne makes 
accretion-induced collapse (AIC) of a WD to a neutron star 
possible, and thus O-Ne WDs, rather than C-O WDs , pref- 



erentially undergo AIC (INomotol 119841: IMivaii et al. 1980; 
Canal et al1ll990l: ICanal et alJll992t iGutierrez et al.l fl99e : 
Gutierrez et al.ll2005l : IWanaio et al J 120091 ) . AIC may result 
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both from accretion following Roche-lobe overflow from a 
non-degener ate companion (the "single degenerate" AIC 
channel; e.g. iNomoto et al.l [19791 : INomoto fc Konddll99ll ). 
as w ell following the m erger of two WDs in a binary system 
(e.g. lYoonet al.l[2u07l ). 

There is no direct observational evidence for AIC and 
so its rate remains uncertain. Numerical simulations of AIC 
find that it is accompanied by a weak explosion (~ 3 x 10 49 
ergs) that synthesizes only a ve ry small amount of N i 
(< 10~ 3 M W in current models') dWooslev fc Baronl Il992l : 
iDessart et al.l 120061 ; lAbdikamalov et al.l 120101 ) , unless there 
is a strong, large-scale magnet i c field in the prog enitor 
white dwarf (|Dessart et al.ll2007l ; iMetzger et al] |2008h The 
neutrino-driven wind from the proto-neutron star does pro- 
duce ~3x 10~ 3 Mq of neutron-rich e lements (specifically , 
matter with an electron fraction < 0.4; IDessart et al.ll2006l ). 
The low total abundance of such elements in the galaxy im- 
plies an AIC rate of <; 1% of the Type la rate, if the calcu- 
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lations of iDessart et al.l (|2006f ) are representative of AIC in 
general. On empirical grounds, lOian fc Wasserburg (|20071 ) 
have suggested that AIC may in fact be the origin of the 
third peak r-process elements. 

Because a WD accretes angular momentum as well as 
mass, it is likely to be rapidly rotating as it approaches M c h- 
Evidence in favor of this expectation is provided by the re- 
cent discovery that the X-ray pulsator RX J0648. 0-4418 is an 
accreting WD wi th a mass > 1.2M^ and a rotational period 
P ~ 13 seconds (jMereghetti et alj :2009). Rapid rotation of 
the WD prior to collapse wou ld make AIC a strong sour ce of 
gravitational wave emission IjAbdikamalov et al,ll2010i ). An 
electromagnetic counterpart to such a detection would sig- 
nificantly increase its scientific yield. More generally, the de- 
tection of AIC would provide important new constraints on 
binary WD evolution and the origin of Type la SNe. New 
optical transient surveys operating with high cadence and 
wide fields of vi e w (e.g., Palomar Transient Factory [PTF; 
iRau et all [20091 ; lLaw et all 120091 ] and the Pan-STARRs-1 
Medium Deep Survey) are sensitive to much fainter tran- 
sients than previous surveys - AIC may thus finally be de- 
tectable. In order to increase the probability of such a de- 
tection, it is critical to have diagnostics that can isolate AIC 
among the deluge of transients that are being discovered. 

In this paper, we present Monte-Carlo calculations of 
the optical-infrared lightcurves and spectra of transients 
produced by AIC. Our radiative tra nsfer calculations ar e 
motivated by the model presented in iMetzger et alj (|2009bh . 
They argued that the observational signature of AIC can, in 
some circumstances, be dominated by the ejecta produced 
by an accretion disk around the newly formed neutron star, 
rather than by the weak SN that accompanies the AIC it- 
self. To mo tivate our radiative tra nsfer calculations, we first 
summarize iMetzger et al.l (]2009bJ)'s calculation of the disk's 
evolution and nucleosynthesis ( qi.ip . We then describe our 
radiative transfer calculations ( i|2.1[) and the composition 
and structure of the disk ejecta ( £|2.2[) . Our main results are 
presented in Sj3] We conclude by discussing the key proper- 
ties of AIC, and comparing them to other predicted faint, 
rapid transients (Q. 

1.1 Summary of the Model 

Depending on the rotation profile of the WD prior to 
colla pse, and in particular the degree of differential rota - 
tion (|Saio fc Nomotdl2004lrYo"on fc Langerll2004l|Pircll200Sft . 
the neutron star formed during AIC is surrounded by 
a rotationally supported disk (|Dessart et al.l 12009 . 120071 ; 
lAbdikamalov et al.ll2010r ). For even modest differential rota- 
tion (that does not significantly change the Chandrasekhar 
mass), the disk mass is Mo ~ 10 -2 — 0.1 M©. However, 
if the WD is uniformly rotating, the disk mass is signif- 
icantly lower, even if th e WD is rotating near break-up 
JAbdikamalov et alJlioiol ). 

Most of the rotationally supported material is initially 
just outside the proto-neutron star at a radius of Rd ~ 100 
km. The disk is born sufficiently hot that nuclei are dissoci- 
ated into free nucleons and weak interaction processes have 
time to change the neutron to proton ratio. Under these con- 
ditions, the composition of the disk can be usefully described 
using the electron fraction 



where n p and n n are the proton and neutron densities, re- 
spectively. Weak interactions in the disk are dominated by 
the pair capture reactions 
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where the resulting neutrino emission is also the domi- 
nant coolant for the disk material because the photons are 
trapped and cannot escape on the timescales of interest. 
Although the disk is formed very hot, electron degeneracy 
suppresses the formation of electron/positron pairs. The re- 
sulting electron excess favors electron captures (eq. [3]) over 
positron capt ures (eq. [2]), thus d riving the composition neu- 
tron rich (e.g. iPruet eti alll2003h . 

As the disk accretes onto the proto-neutron star, the 
remaining mass spreads to larger radii in order to con- 
serve angular momentum. As the disk spreads, the tem- 
perature drops (for the majority of the mass) and thus 
the neutrino cooling and weak reactions become less ef- 
ficient. When the bulk of the disk reaches Rd ~ 1000 
km, the disk has cooled to T ~ 10 10 K , and weak inter- 
actions freeze out l|Metzger et al.ll2009al ). Absent a central 
neutron star - as is the case for analogous disks formed 
during binary neutron star mergers - weak interactions 
freeze out when the material is still quite neutron rich, with 
Y e ~ 0.2 — 0.4. During AIC, however, the proto-neutron 
star emits copious neutrinos in the fi rst few seconds after 
formation (|Burrows fc Lattimerlll986T ). which irradiate the 
disk as it spreads. Because the neutron star is deleptoniz- 
ing, electron neutrinos dominate the emission and the rate 
of neutrino captures v E + n — > p + e~ becomes comparable 
to the rate of pair captures (eqs. [2]-[3]). For a significant 
fraction of the mass, this increases the electron fraction at 
freeze-out to Y e > 0.5. 

As the disk continues to spread, the free nucleons in the 
disk recombine to form heavy nuclei once T ~ few xlO 9 K. 
The formation of He alone releases ~ 7 MeV per nucleon, 
which is larger than the gravitational binding energy of the 
disk when fusion commences (~ 2 MeV per nucleon). As a 
result, the disk is unbound, with its slow viscous spreading 
transformed into a supersonic outflow that reaches u~0.1c. 
The end state of the fusion that initiates the disk's destruc- 
tion is determined by the density, temperature, and electron 
fraction (the latter because weak interactions have frozen 
out); the parameters are such that the ejecta is nearly in nu- 
clear statistical equilibrium (NSE). Since there are roughly 
equal numbers of protons and neutrons, a significant frac- 
tion of the unbound disk fuses to form 56 Ni, the decay of 
which can power a supernova- like transient once the ejecta 
expands sufficiently that photons may diffuse out. In the 
remainder of this paper, we calculate the lightcurves and 
spectra of this transient. 

2 ASSUMPTIONS & MODEL PARAMETERS 

2.1 Radiative Transfer Calculations 

We calculate lightcurves and spectra using the time- 
dependent radiative transfer code SEDONA, a detaile d de- 
scription of which can be found in iKasen et al.l (|2006l ) . SE- 
DONA solves the radiative transfer problem in a homolo- 
gously expanding atmosphere described by a velocity-time 
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grid. We include the gamma ray energy emitted from the 
decay chain 56 Ni — > 56 Co — > 56 Fe and evolve the composi- 
tion accordingly. We modified SEDONA to treat 56 Ni as dis- 
tinct, allowing us to include other Ni isotopes in the opacity 
without evolving their composition or including th e ir dec ay 
energy. All line information is from lKurucz fc Belli (|1995T ). 

SEDONA tracks the emitted gamma rays from radioac- 
tive decay and determines the resulting heating through 
Compton scattering and photoelectric absorption. The opac- 
ities are first calculated using the input density and compo- 
sition, the calculated energy deposition, and an initial tem- 
perature guess. The dominant opacities are electron scatter- 
ing and bound-bound line transitions. Our calculations as- 
sume local thermodynamic equilibrium (LTE) , although SE- 
DONA supports non-LTE opacity calculations as well. The 
source function for the lines is calculated in the two-level 
atom formalism. To account for line scattering, the source 
function, 5a, deviates from its LTE form of the Planck func- 
tion, B\, through 



Table 1. Mass Fractions in the Three Fiducial Models. 



Sx = (l-e)Jx+eB x 



(4) 



where the degree of line scattering is encapsulated in the 
thermalization parameter, e, which equals the ratio of ab- 
sorption opacity to tot al opacity. We fixed e = 0.3, motivated 
by the calculations of iKasen et al.l (|2006r ). Given the opac- 
ity and ionization states in LTE, the temperature profile is 
determined from the condition of radiative equilibrium, bal- 
ancing the rate of thermal emission with the rate of photon 
absorption and radioactive energy deposition. The synthetic 
spectrum is determined using a Monte Carlo calculation. 
The spectrum and temperature profile are iterated on un- 
til the temperature structure converges. The final synthetic 
lightcurves and spectra are determined on the last iteration 
by recording the photons that escape the outflow. We set 
the velocity and time bin resolutions in our calculations to 
200 km s _1 and 0.1 day, respectively. 



2.2 Outflow Properties 

Although the outflows produced during AIC are likely to 
be asymmetric (albeit axisymmetric) , we use a ID spheri- 
cally symmetric model to study the lightcurves and spectra. 
Given current uncertainties in the properties of the unbound 
accretion disk and the interaction between the disk and the 
AIC supernova, we believe that this approximation is a rea- 
sonable first step. 

The bulk speed of the outflow is ~ 0.1 c, determined by 
the nuclear energy released as the disk is unbound. Given the 
likely order unity fluctuations about this fiducial speed, we 
set the density distribution of the homologously expanding 
outflow to be 



P = 



ae v/v " if 0<v< 0.05 c 
bv~ c if 0.05 c sC v < 0.2 c 



(5) 



where C parameterizes the density profile of the outflow. 
The constants a and b, which in general are functions of 
time, are determined by continuity of density at v = 0.05 c 
and by fixing the total outflow mass, Mtot- The e- folding 
velocity, « e , is included to ensure that there is not an arti- 
ficially abrupt cutoff in the density at small velocity /radii 
(which can affect the photosphere at late times), but the 











Y e 


range 






Element 


0.425-0.55 


0.45-0.55 


0.5-0.55 




Ca 


3.60 


X 


10~ 5 


4.55 


X 


10~ 5 


1.68 x 10~ 


-5 


Sc 


1.58 


X 


IO" 6 


1.99 


X 


IO" 6 


3.62 x 10~ 


-6 


Ti 


5.63 


X 


10" 4 


2.55 


X 


IO" 5 


4.44 X 10~ 


-5 


V 


1.40 


X 


10~ 6 


1.77 


X 


HT 6 


3.22 x IO" 


-6 


Cr 


2.24 


X 


IO" 3 


1.16 


X 


10" 4 


7.61 x 10- 


-5 


Fe 


1.87 


X 


io- 2 


1.12 


X 


IO" 2 


1.07 x 10- 


-2 


Co 


2.95 


X 


10~ 3 


3.73 


X 


10~ 3 


1.76 x 10" 


-3 


Ni 


8.59 


X 


HT 1 


9.00 


X 


io- 1 


9.87 x 10- 


-1 


Cu 


5.03 


X 


10~ 3 


4.52 


X 


10- 3 


2.21 x 10- 


-4 


Zn 


1.03 


X 


io- 1 


7.31 


X 


10~ 2 


5.90 x 10- 


-5 


Ga 


1.42 


X 


HT 4 


1.40 


X 


IO" 4 


0.00 




Ge 


3.86 


X 


IO" 3 


4.23 


X 


io- 3 


0.00 




As 


5.60 


X 


IO" 6 


7.08 


X 


IO" 6 


0.00 




Se 


5.97 


X 


10" 4 


5.42 


X 


10~ 4 


0.00 




Rb 


1.69 


X 


HT 4 


1.25 


X 


10- 5 


0.00 




Sr 


2.86 


X 


KT 3 


1.12 


X 


10- 3 


0.00 




Y 


1.84 


X 


10" 4 


2.19 


X 


10" 4 


0.00 




Zr 


9.46 


X 


10" 4 


1.20 


X 


IO" 3 


0.00 




Mo 


3.40 


X 


10~ 6 


4.30 


X 


IO" 6 


0.00 





Table 2. Mass Fraction of 56 Ni in the Three Fiducial Models. 



0.425-0.55 



Y e range 
0.45-0.55 



0.5-0.55 



5.38 x 10" 1 6.49 x 10~ 



9.97 x 10- 1 



exact value is not crucial; we take v e — 500 kms -1 . Equa- 
tion © implies that the bulk of the mass in the outflow has 
0.05c ^ v < 0.2c. The precise value of the density power 
law index C within this velocity range does not significantly 
affect the observational signatures of AIC. We will focus on 
C = 4 but we carried out calculations for C — 2 — 4, which 
covers the full range of possibilities, from most of the mass 
being at high velocity (C = 2) to most of the mass being 
at smaller velocity (C = 4). The lightcurves in the models 
with smaller C peak somewhat earlier (At = 0.2 d between 
C = 2 and C = 4) and fall more rapidly after the peak, 
as would be expected since more of the mass is at high ve- 
locity. Nonetheless, the peak bolometric luminosity remains 
nearly the same. The spectra at early times are relatively 
independent of C, but there are some differences at t ;> 3 
days when individual spectral features begin to become more 
prominent. We will highlight these differences below, but we 
focus on the case of C = 4. 

Metz ger et al.l (|2009br ) showed that different models for 
the evolution of the accretion disk around the NS (e.g., dif- 
ferent values of the Shakura-Sunyaev 1973 viscosity param- 
eter a) can lead to somewhat different composition in the 
resulting ejecta. Specifically, the electron fraction is within 
the range 0.4 < Y e < 0.65, but the distribution of mass as 
a function of Y e depends on the details of the disk's evolu- 
tion. To cover this parameter space we studied three fidu- 
cial models which include an increasingly large amount of 
neutron-rich material at low F e : (1) Y e — 0.5 — 0.55, (2) 
Y e = 0.45 - 0.55 and (3) Y e = 0.425 - 0.55. The models 
all have equal mass per unit Y e . We believe that the mod- 
els with the larger range of electron fraction (in particular, 
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Y e = 0.42 5 — 0.55) are the mos t realistic because the disk 
models in lMetzger et al.l l|2009bl ) in general have a range of 
Ye at different radii and times. However, the differences be- 
tween these three different models give an indication of the 
systematic uncertainties in our predictions introduced by 
uncertainties in the outflow composition. Lacking a hydro- 
dynamic calculation of the disk evolution and its transition 
to the ejecta phase, we take the composition of the ejecta to 
be uniform at all radii/velocites. 

Given the electron fraction in the ejecta, we deter- 
mined the composition of each model using the n ucle- 
osynth esis calculations of IWooslev fc Hoffman! fl992j) and 
Scitcnzahl et al. (2008), averaging over the range of Y e in our 
model , and using the solar system abundances from I Arnettl 
(1996) to determine absolute abundances. Table [T] summa- 
rizes the composition of each of the models and Table [2] 
gives the mass fraction of 56 Ni in each model. In determin- 
in g the composition of our mo dels, we used the calcu l ations 
of IWooslev fc Hoffman! l|l992J ) and ISeitenzahl et~aT1 (|2008l ) 
with thermodynamic properties most similar to AIC ejecta. 
Specifically, for 0.425 ^ Y e < 0.5, we used the finite entropy 
abundances calculated from the a-pro cess with T = W 10 K 
and p = 2 x 10 7 g cm~ 3 in Table 3 of IWooslev fc Hoffman! 
(1992), corresponding to a mild a-rich freeze out. The en- 
tropy and expansion time scale used in these calculations 
are similar to those expected in the disk when a— particles 
form and outflows begin. For 0.5 < Y e ^ 0.55, we use the 
abundances calculated from NSE with T = 3.5 x 10 9 K an d 
p = 10 7 g cm" 3 in Figure 3 of ISeitenzahl et all J2008I ). 
Because AIC outflow s likely possess a m oderate entropy 
~ 3— lOfcs baryon -1 l|Metzger et al.ll2009bl ). and hence may 
u ndergo a mild a— rich fre eze-out, the pure NSE calculations 
in ISeitenzahl et a.1.1 |2008) do not stricly apply. However, by 
ex trapolating the res ults of full nucleosynthesis calculations 
bv lPruet et al.l (J200J) to the relevant entropy regime, we find 
that the yield of 56 Ni, the dominant isotope in this compo- 
sition range, is well-approximated by its NSE abundance. 

Since the nucl e i are in their post-decay phase in 
IWooslev fc Hoffman! (1 19921 ). we converted the entire frac- 
tion of 56 Fe to its parent 56 Ni. Although some matter with 
Y e < 0.425 may be present in AIC ejecta, we did not in- 
clude it in our models due to the u ncertain and con t radic - 
tor y nucleosynthesis i nform ation in ISeitenzahl et al.l ([2008) 
and lHartmann et al.l (| 19851 ). We also did not include Br and 
Kr in our calculatio ns due to lack of line information in 
iKurucz fc Belli (1 19951 ). Finally, our calculations only include 
heating due to the decay of 56 Ni, although other radioac- 
tive isotopes are produced and may contribute to the light 
curve. In particular, we find that 57 Ni (half-life T\n = 35.6 
hr) and 62 Zn (ti/2 = 9.2 hr) produce specific heating rates 
that are a factor ~ 6 and ~ 12 times higher, respectively, 
than 56 Ni on timescales t <C T\i%. However, in most cases the 
much lower mass fractions of these isotopes relative to 56 Ni 
implies that their contribution to the total heating is rela- 
tively minor. One exception is the mass fraction of 62 Zn in 
our Ye = 0.425 — 0.55 model, which we estimate is only a fac- 
tor ~ 4 less than 56 Ni and thus could contribute appreciably 
to the very early light curve. 

As discussed in 91.11 the total mass in the ejecta 
produced during AIC is uncertain, and d epends on the 
differen tial rotation in the WD progenitor. iMetzger et al.l 
l|2009bl ) determined the mass of unbound disk material to 



be a few tim e s 10~ Mq using collapse calculations from 
iDessart et al.l (J2006J, 120071 '). but the ejecta mass could be 
significantly lower if the WD is more uniformly rotating 
(Abdikam alov et al.l 120101 ). We studied a range of total 
masses for the ejecta but focus on Mtot =3x 1O _3 M0 and 
Mtot — 10 _2 Mq here because this is the range of masses 
for which the disk ejecta produces an observationally inter- 
esting luminosity. The lightcurves and spectra have similar 
behaviour for different total masses, simply shifted to some- 
what earlier times and shorter wavelengths for lower total 
mass, as would be expected u sing analytic estimates based 
on Arnett's law (|Arnettlll98"2l ). 

3 RESULTS 

Figure [T] shows the bolometric and R-band lightcurves pro- 
duced by AIC ejecta for the three different compositions 
summarized in Tables [T] and [21 and for two different to- 
tal masses. The lightcurve evolution is consistent with an- 
alytic expectations b ased on the c ompetition between ex- 
pansion and diffusion |Arnettll 19821 ) . The low mass and high 
speed of the ejecta imply that the rise time is extremely 
fast ~ day; the subsequent decay is somewhat slower, with 
the luminosity decaying significantly after ~ 4 — 5 days for 
Mtot = 10" 2 M Q and ~ 2 - 3 days for Mtot = 3 x 10" 3 M Q . 



The peak bolometric luminosity is 



ergss 



for Mtot — 10~ 2 Mq and a factor of ~ 3 lower for Mtot =■ 
3 x 10~ 3 Mq. This is a factor of ~ 100 fainter than a typ- 
ical Type la supernova, which is consistent with the much 
smaller Ni masses produced during AIC. 

For a given total mass, Figure [T] shows that the models 
with smallest range in electron fraction (Y e = 0.5 — 0.55) are 
the brightest and broadest. The dependence on the compo- 
sition is, however, modest and is reasonably degenerate with 
a different total mass for the ejecta. This implies that the 
uncertainty in the detectability of AIC is due primarily to 
uncertainty in the total mass of the disk ejecta, not the com- 
position of the material. The dependence on Y e in Figure [T] 
arises because only matter with Y e > 0.48 produces signif- 
icant Ni, which both powers the lightcurve and contributes 
to the opacity of the ejecta. Lower Y e material ends up as 
higher Z elements which do not contribute to the radioactive 
heating (note that for even lower Y e ~ 0.1 — 0.3, the decay 
of neutron -rich ejecta can prod uce a heating rate similar to 
that of Ni; IMetzger et"alll2010| ). 

The bumps in the AIC lightcurves in Figure [T] have the 
same physical cause as thos e seen in som e Type la SNe (see 
IPinto fc Eastman! |2000| and lKasenll2006l for a discussion of 
the latter). As the ejecta expands, the matter cools, and 
ions begin to recombine. Singly ionized ions have a greater 
number of lines in the optical and infrared part of the spec- 
trum, and are therefore more efficient at fluorescing photons 
from bluer to redder wavelengths. Because the ejecta is more 
transparent at redder wavelengths, this redistribution effect 
allows the radiation energy stored in the ejecta to escape 
more rapidly, leading to a rise in the bolometric light curve. 

Figure [2] shows the U-B, B-V, B-J, and R-I colors as 
a function of time for the Mtot — 10 _2 Mq models from 
Figure [l] We define the colors here using the ratio of the 
luminosities \L\ at the centers of the appropriate bands, 
with Au = 3611 A, A B = 4611 A, A v = 5851 A, Ar = 
7011 A, Ai = 8111 A, and Aj = 12351 A. Figure [2] shows 
that AIC is predicted to become increasingly red in the first 
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Figure 1. Predicted bolometric and R-band lightcurves of AIC transients, calculated for different values of the total ejected mass Mtot- 
The R-band center was taken as Ar = 7011 A. Thick and thin lines correspond to Mtot = 10 _2 Mq and 3 X 10 — 3 Mq, respectively. 
Different linestyles correspond to different ejecta compositions, specifically different assumptions about the range of electron fraction in 
the ejecta (see Tables [TJ and [2]l . 



~ day, after which there is not as much color evolution as a 
function of time. This is particularly true for our preferred 
models which have a wider range of electron fraction in the 
outflow. The model with the narrow range of Y e (0.5-0.55) 
is bluer and shows somewhat more color evolution at early 
times, as a result of being almost entirely 56 Ni (Tables[l]and 
0). Note, however, that although the emission as a whole be- 
comes redder with time, the optical colors are slightly blue, 
with B - V ~ -1 to -2. 

The lower mass models with Mtot = 3 x 10~ Mq show 
the same qualitative U-B color evolution as the M to t = 
1O _2 M0 case. In particular, they reach the same color 
plateau, but somewhat faster than in the high mass case 
because the ejecta becomes transparent earlier. One differ- 
ence, however, is that for Mtot = 3 x 1O _3 M0, the B-V and 
B-J colors become redder at day 3, and the R-I colors peak 
at the same value 2 days earlier than the high mass cases and 
rise and fall faster. Models with shallower density profiles (C 
= 3 and 2 in eq. [5]) show the same qualitative behavior for 
U-B, B-V, and U-J as our standard C = 4 model; they are, 
however, slightly bluer in U-B throughout (by ~ 0.6, 1.2), re- 
spectively, while in B-V and B-J they are equal until day 3, 
after which they overturn and become redder by ~ 1,2 mag- 
nitude, respectively. In R-I, they peak at the same value but 
~ 0.9, 1.5 days earlier, respectively, and rise and fall faster. 

Our predicted color evolution may be sensitive to non- 
LTE effects. In particular, the red colors are primarily due 
to the fact that in these models a significant fraction of the 
mass is in iron and nickel even at the highest velocities. The 
outer layers of the ejecta are cooler, and in a lower ioniza- 
tion state, and so they very strongly blanket the bluer wave- 
lengths. Non-thermal ionization from the decay products of 
56 Ni could in principle keep the outer layers of the ejecta 
more ionized. This would reduce the degree of line blanket- 
ing in the blue part of the spectrum (and the fluorescence to 
redder wavelengths) leading to an overall less red spectrum. 
This possibility will be studied in future work. 



The spectra at 1 day and 5 days are shown in Fig- 
ure[3]for the Mtot — 1O~ 2 M0 models. Near the peak of the 
lightcurve (~ 1—3 days), the emission peaks at ~ 3000—6000 
A and is reasonably red, consistent with Figure [3 At 1 day, 
individual spectral features are not discernible: given the 
high speeds, the lines have Doppler broadened and merged. 
Spectral lines do become more apparent at later times. In 
particular, at ~ 3 — 5 days the calcium triplet at ~ 9000 
A produces a luminosity comparable to the total luminos- 
ity of the ejecta, despite the fact that the mass fraction of 
Ca is ~ 3 x 10~ 5 . This is the most prominent individual 
spectral feature seen in any of our calculations, although we 
note that the actual strength of the Ca II line may be some- 
what overestimated due to our LTE assumption. The spectra 
are generally Ni dominated but Doppler broadening makes 
it very difficult to identify individual spectral features. Al- 
though the results in Figure [3] are for M to t — 1O _2 M0 and 
p oc v~ 4 in equation (0), the general conclusions about the 
spectra of AIC drawn here are reasonably robust and also 
apply to the other calculations we have carried out. In par- 
ticular, the lower mass models with M to t — 3 x 10 _3 Mq 
have similar spectra at day 1. Models with different density 
profiles have similar spectra at day 1, but steeper velocity 
profiles (larger C) retain the Ca II triplet line and the peak 
at ~ 5000 A to later times, simply because more mass is 
concentrated at lower velocities. 

To assess the potential impact of the weak SN accom- 
panying AIC on our predicted lightcurves and spectra, we 
also carried out a calculation in which we included the likely 
nucleosynthetic products of such a SN in our radiative trans- 
fer calculation. We assume that the SN material is uniformly 
mixed with the disk wind ejecta, which ignores the poten- 
tially interesting dynamics of the interaction between the 
disk wind and the SN. We set the mass in the SN ejecta to 
be 3 x 10~ 3 M ( 7), mo tivated by hydrodynamic simulations in 
iDessart et al.l ([20061 ) - note that since there is no significant 
outer stellar envelope in AIC, the "SN" mass is largely the 
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Figure 2. U-B, B-V, B-J, and R-I color evolution for the model with ejecta mass Mtot = 10 — 2 Mq. The models with the wider range 
of electron fraction are, we believe, more physical (Y e = 0.425 — 0.55 and Y e = 0.45 — 0.55); see text. These initially become redder with 
time in U-B, with little evolution after a few days. The spectra at 1 and 5 days are shown in Figure[3] The colors arc defined using the 
ratio of the luminosities \L\ at the following wavelengths: Atj = 3611 A, Ab = 4611 A, Av = 5851 A, Ar = 7011 A, Ai = 8111 A, and 
Aj = 12351 A (chosen to match the wavelength bins used in the radiative transfer calculation). 



mass unbound by the neutrino-driven wind from the proto- 
neutron star in the first ~ 1 sec after collapse. Given tha t 
most of this ejecta has Y e ~ 0.25 - 0.4 (JDessart et al.ll200ej ). 
it will primarily form elements near the second r— process 
peak at A ~ 130 JHoffman et al.lll997t ). Although a detailed 
nucleosyntheis calculation would be necessary to determine 
the wind's precise nuclear yield, for simplicity we include in 
our calculation elements in the range of the second peak (Z 
= 45-60) with an abundance distri bution norma lized to the 
observed solar values (taken from lArnettlll99q ). We omit- 
ted I and Xe becaus e no line information is available in 
iKurucz fc Belli (|1995l ). The resulting model has a slightly 
larger peak luminosity and rise time as compared to our 
fiducial model without the additional matter (shown in Fig. 
[TJ , but the lightcurve is somewhat broader and has a slower 
decay. This is simply a consequence of the larger total mass 
of the ejecta. Although the AIC plus SN wind case may be 
somewhat redder than AIC alone at late times, we did not 



find any signature (e.g., spectroscopic) that could distin- 
guish entrained SN material from a disk wind with a some- 
what larger mass. We note, however, that the available line 
data on very heavy elements is incomplete, which may par- 
tially explain the lack of discernable spectroscopic features. 



4 DISCUSSION 

Although the theoretical possibility of AIC has been recog- 
nized for decades (| Canal fc Schatzman 1976:Nomot o et al.l 
Il979f ). only recently have transient surveys had sufficient 
depth, field of view, and cadence to potentially detect 
and characterize these events. The identification of AIC, 
or the ability to place meaningful constraints on its rate, 
would have important implications for a variety of astro- 
physical problems, including degenerate binary evolution, 
the progenitors of Type la SNe, and the fo rmation chan- 
nels of globular cluster NSs (|Grindlavlll987l ). low mass X- 
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Figure 3. Spectra for the model with ejecta mass Mtot = 10 — 2 Mq. At early times the spectra are dominated by Doppler-broadened 
Nickel features, with no individual spectral features identifiable. At 5 days, however, the Ca triplet at 9000 A is quite prominent. This 
line is generically present at late times (~ 2 — 5 days) in all of our models. 



ray binaries (Ivan den Heuve]||l984T ). and millisecond pulsars 
(|Grindlav fc Bailvnlll98St ). 

Because AIC transients are powered by the decay of 
56 Ni, they have many qualitative similarities to Type la SNe. 
Although the observed lightcurves are too dim to perform 
cosmology, characterizing the correlatio n between pe ak lu- 
minosity and decay timescale (i.e., the IPhillipsI Il993l rela- 
tion) for AIC would be useful for observationally identify- 
ing these events. For example, the duration and peak lu- 
minosity of AIC transients increases with the ejecta mass, 
which probably incre ases in proportion to the initial mass 
of the accretion disk (jMetzger et al.ll2009bl ). The disk mass, 
in turn, is set by the angular mome ntum and rotational 
profile of the WD prior to collapse (Des sart et al.l I2006J ; 
lAbdikamalov et al.l 120101 ) , which also determines whether 
AIC is expected to be a powerful source of gravitational 
waves (|Otg|2009j). Everything else being equal, more lumi- 
nous AIC transients are thus more likely to be accompanied 
by a detectable gravitational waves, making AIC transients 
a promising electromagnetic counterpart to events detected 
with ground-based interferometers such as Advanced LIGO. 

We now briefly address the prospects for the "blind" 
detection of AIC with current and future transient surveys. 
Although the rate of AIC is uncertain, the amount of rare 
elements produced by the neutron-rich SN ejecta constrains 
the rate 7?aic to be < 10~ 4 yr -1 per Milky Way-type 
galaxy (or <; 10~ 2 of the S N la rate) if the AIC calcula- 
tions of lDessart et al.l (|2006r ) are representative. Assuming 



that th e AIC rate is proportio nal to the blue stellar lumi- 
nosity (jKopparapu et al.l 120081 ) ■ a Galactic rate of i?Aic = 
10 -4 .R_4 yr -1 corresponds to a volumetric rate of 10 _6 .R_4 
Mpc -3 yr -1 . For an ejecta mass of M e j = 10 _2 Mq, we find 
that AIC results in an optical transient with a peak lumi- 
nosity Lpcak = few xlO 41 ergs s" 1 (Fig. [T]). The resulting 

1/2 

luminosity distance is 600(400) [100]L 4 i Mpc for a magni- 
tude limit of 25(24)[21], where L41 = L pcak /10 41 _ergs_sJ 



The PanSTARRs Medium Deep Survey (MDS; lKaiser et all 
[2002J), which covers ~ 84 deg 2 in g and r down to AB mag- 



nitude ~ 25, should thus detect ~ 2L 41 R-4 AIC transients 
per yr. The 5-d ay cadence surve y of the Palomar Transient 
Factory (PTF; lLaw et al.l 120091 ) . which surveys an active 
area ~ 2700 deg 2 to a limiting AB magnitude of 21, will 
detect 0.3Ll{ 2 R- 4 yr" 1 . Thus, for i? A ic ~ 10" 4 yr -1 and 
L41 ~ few we conclude that MDS and PTF should detect 
~ 1 — 3 AIC per year. As we discuss below, identifying 
AIC transients for follow-up observations will be challeng- 
ing, so these estimates are probably optimistic and motivate 
faster c adence surveys such as the Palomar P60-FasTING 
survey ( JKasliwal et al. 2010). Prospects for detection are 
better with the La rge Synoptic Survey Telescope (LSST; 
IStrauss et al.l I2010T ) . which will image the entire sky down 
to a limiting magnitude ~ 24.5 every 3—4 nights and should 
detect AIC events at a rate ~ 10 3 L 4 ( i?_4 yr -1 . Unfor- 
tunately, the fraction of AIC that have ejecta masses of 
~ lO~ 2 Af0, and thus peak luminosities of L pca k = few xlO 41 
ergs s , is very difficult to predict given uncertainties in 
the differential rotation of the progenitor WD. Absent addi- 
tional theoretical progress, it will thus not be easy to turn 
upper limits from transient surveys into constraints on the 
rate of AIC. 

Although the prospects that upcoming transient sur- 
veys will detect AIC are promising, their identification will 
be challenging due to the deluge of other transients that 
will be discovered (and, indeed, are already being discov- 
ered). Other thermal transients predicted to occur on ~ 
few d ay timescales inclu de the failed deflagrations of C/O 
WDs |Livne et al.ll2005l ), ".la" SNe due to unstable ther- 
monuclear He flashes from WD binaries (|Bildsten et al.l 
120071 ; IShen et al.ll2010|}. shock breakout from Type la su- 
pernova ( Piro et alj|2010 j ). and radioactively -powered tran- 
sients from NS mergers (jMetzger et al.ll2010T ). 

Although NS mergers and other WD-related explosions 
may originate from a similar stellar population as AIC, a 
defining characteristic of AIC is that the majority of the 
ejected material is processed through NSE. This means that 
few low or intermediate-mass elements are likely to be pro- 
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duced (such as He, O, C, and Ti), in contrast to what is 
expected from .la's or the failed detonation of a WD. We 
note, however, that trace amounts of Ca and Ti are pro- 
duced in some of our models (see Table [T]) , which still may 
produce strong lines despite their low abundances (e.g. the 
Ca triplet at ~ 9000 A; Fig. [5J. A perhaps more dramatic 
difference between AIC transients produced by a disk around 
the newly formed NS and other WD-related events are the 
large outflow speeds ~ 0.1c in AIC, which leads to highly 
velocity-broadened lines. 

The speed of the ejecta from NS mergers, as in AIC, 
is expected to be large. A distinguishing feature in this 
case, however, is that NS merger ejecta is more neutron-rich 
than AIC ejecta, because in AIC the accretion disk material 
that ultimately becomes unbound is first irra diated by elec- 
tron n eutrinos from the proto-neutron star (JMetzger et al.l 
l2009al lbr). Although the current line data for the r— process 
nuclei produced in NS merger ejecta is insufficient to make 
concrete spectroscopic predictions, it is safe to say that NS 
mergers are unlikely to resemble any spectral class of SNe yet 
detected. AIC and NS merger transients may also be distin- 
guished by their light curve evolution. In Figure [4] we com- 
pare the bolometri c light curves from A IC (this paper) and 
NS mergers ffrom iMetzger et al.l 120101 ) assuming an ejecta 
mass M e j = 10 _2 Mq in both cases. Although the overall 
brightness and late-time evolution is similar, the NS merger 
light curve shows an earlier, more pronounced peak. This 
difference results from the different radioactive heating rate 
in the two cases. AIC transients are powered by the decay 
of 56 Ni, which produces an approximately constant heating 
rate prior to the 56 Ni half-life ~ 6 days. NS merger tran- 
sients, by contrast, are powered by a large number of de- 
caying nuclei with a wide distribution of half-lives, resulting 
in a total heating rate that decreases approximately as a 
power-law E oc i -1 ' 3 . NS merger transients peak ear- 
lier simply because the radioactive heating is concentrated 
at early times. Distinguishing between NS mergers and AIC 
based on light curve evolution alone would thus probably 
require catching the event prior to maximum. 

Several sub-luminous, rapidly-evolving Type I su- 
pernovae have recently been discovered, which could in 
principle be r elated to AIC. These i nclude SN 2008h a 
(IValenti et al.l 120091; [Foley et all [20091 ; iFolev et all l2010h 
2005E (|Perets et al.l [20091 ). and 2002b i l|Poznanski et all 
120101 ). In all of these events, the low 56 Ni masses inferred 
from their peak luminosities (~ 10 -3 — O.IMq) are largely 
consistent with the quantity predicted from AIC. However, 
the presence of low and intermediate mass elements (e.g. Ca, 
C, O, He) in the spectra of these events with relatively low 
velocities (~ 2000 — 10000 km s _1 ) rules out the standard 
Ni-rich AIC scenario outlined in thi s paper . 

As discussed by IMetzger et al.l (|2009bl ). the timescale, 
velocity, and composition of the AIC transient can be modi- 
fied if a substantial quantity of WD material remains at large 
radii following AIC, either because the WD is still collapsing 
or because some matter remains centrifugally supported in a 
remnant disk cre ated by a super-Mch WD- WD merger (e.g. 
lYoon et al.ll2007l ). In the latter case, depending on the mass 
and composition of the WDs, up to ~ O.IAIq in C, O, Ne, 
or He could remain at a radius ~ 10 9 cm. Once the SN ex- 
plosion or Ni disk wind impacts the remnant material, it will 
(1) shock heat the material to a few times 10 9 K, which may 
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Figure 4. Comparison of the bolometric light cur ves of the tran- 
sient s from AIC (this paper) and NS mergers IIMetzger et al.l 
[2010J). Both models were calculated for a total ejecta mass of 
Aftot = W~ 2 Mq. The NS-NS merger transient is produced by 
the decay of a large number of neutr on-rich nuclei created by 
r-process nucleosynthesis in the ejecta; IMetzger et al.l 11201011 es- 
timated a thermalization efficiency of etherm ~ 0-5 f° r these de- 
cays. This figure highlights that distinguishing NS merger and 
AIC transients using lightcurve evolution alone will likely require 
early-time observations at < 1 day. 



synthesize intermediate mass elements but will not disasso- 
ciate the Ni and will leave some unburned WD material; (2) 
slow the ejecta from v ~ 0.1 c to a few thousand km s _1 . 
A combination of slower ejecta and higher opacity would 
lengthen the duration of the light curve, making such an 
event easier to detect. Such an "enshrouded AIC" scenario 
may therefore still be capable of explaining events such as 
SN 2008ha, 2005E, and 2002bj. In order to make a quanti- 
tative comparison to observations, however, hydrodynamic 
simulations will be required to determine the dynamics and 
nucleosynthetic yield of the complex interaction between the 
prompt SN, the Ni-rich disk wind, and the remnant WD ma- 
terial. Moreover, the timescale of AIC relative to the WD- 
WD merger itself is uncertain; any significant delay would 
likely lead to less debris at large radii and a "cleaner" event 
analogous to that calculated in this paper. 

On this note it is also important to distinguish between 
AIC due to accretion from a binary companion, as consid- 
ered in this paper, and electron-capture supe rnova in AGB 
stars with O-Ne-Mg cores (for recent work, see lKitaura et al.l 
2006; Wanaio et alj2009f ). Although the collapse mechanism 
in these two cases is similar, an important difference in the 
AGB case is the presence of an extended hydrogen enve- 
lope. This additional overlying material both increases the 
56 Ni yield of the explosion and ext ends the duration o f the 
supernova due to its larger opacity (IWanaio et alj|2009f ) . Al- 
though the lack of hydrogen in the spectra of SNe 2008ha, 
2005E, and 2002bj rules out electron capture SNe in these 
events, electron capture events may ex plain slower evolv- 
ing sub-luminous SN such as 2008S (e.g. iPrieto et aI1l2008l ; 
iThompson et al.ll2009l ; IPumo et al.ll2009h 

Our calculations could be improved in a number of 
ways. In particular, future work will require more realistic 
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multi-dimensional models in which the composition is calcu- 
lated from a hydrodynamic simulation of the disk evolution 
beginning with free nucleons. This is important for deter- 
mining the final composition as a function of radius, which 
in turn can influence the spectroscopic predictions. 
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